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Abstract 

We consider the dynamics of a periodic chain of N coupled overdamped particles un- 
' der the influence of noise. Each particle is subjected to a bistable local potential, to 

a linear coupling with its nearest neighbours, and to an independent source of white 
I noise. The system shows a metastable behaviour, which is characterised by the loca- 

p I ■ tion and stability of its equilibrium points. We show that as the coupling strength 

increases, the number of equilibrium points decreases from 3^ to 3. While for weak 
coupling, the system behaves like an Ising model with spin-flip dynamics, for strong 



^ ' coupling (of the order iV^), it synchronises, in the sense that all particles assume al- 

most the same position in their respective local potential most of the time. We derive 
the exponential asymptotics for the transition times, and describe the most probable 
^Sj ■ transition paths between synchronised states, in particular for coupling intensities be- 

^ I low the synchronisation threshold. Our techniques involve a centre-manifold analysis 

. of the desynchronisation bifurcation, with a precise control of the stability of bifur- 

cating solutions, allowing us to give a detailed description of the system's potential 
landscape. 
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H ■ 1 Introduction 

, 

Lattices of interacting deterministic multistable systems display a wide range of interesting 
behaviours, due to the competition between local dynamics and coupling between different 
sites. While for weak coupling, they often exhibit spatial chaos (independent dynamics 
at the different sites) , for strong coupling they tend to display an organised collective be- 
haviour, such as synchronisation (see, for instance |BM96l [CMPVV96t [JohOTj INMKV97] . 
and [PRKOli ICFOS] for reviews). 

An important problem is to understand the effect of noise on such systems. Noise 
can be used to model the effect of unresolved degrees of freedom, for instance the in- 
fluence of external heat reservoirs (see, e.g., |FKM65l ISLTTl IEPRB99[ IRBTOOl IRBT02j ). 
which can induce currents through the chain. The long-time behaviour of the system is 
described by its invariant measure (assuming such a measure exists); however, for weak 
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noise, the dynamics often displays metastability, meaning that the relaxation time towards 
the invariant measure is extremely long. 

Metastability has been studied extensively for particle systems with stochastic dynam- 
ics. In these models, the transition from one metastable state to another usually involves 
the gradual creation of a critical droplet through small fluctuations, followed by a rapid 
transition to the new state. The distributions of transition times, as well as the shapes 
of critical droplets, have been investigated in detail (see in particular |dH041 lUVOSj for 
reviews, and references therein). 

In these lattice models, the local variables take only a finite number of discrete values, 
which are independent of the interaction with other sites. In the present paper, we consider 
by contrast a model with continuous on-site variables. This leads to a system of interacting 
stochastic differential equations (also called interacting diffusions, see, for instance, |DG88j 
for a review of asymptotic properties in the mean- field case). It turns out that while this 
system has a similar behaviour to stochastic lattice models for weak coupling, the dynamics 
is totally different for strong coupling: There are only 3 equilibrium configurations left, 
while the activation energy becomes extensive in the number N of particles. For large A^, 
the system's behaviour is closer to the behaviour of a Ginzburg-Landau partial differential 
equation with noise (see, e.g. [EII01l[Rou02] ). The transition between the strong-coupling 
and the weak-coupling regimes involves, as we shall see, a sequence of symmetry-breaking 
bifurcations. Such bifurcations have been studied, for instance, in |QC04| for the weak- 
coupling regime, and in |McN99[ IMcNn2[ IWat93a[ IWat93bj for systems of coupled phase 
oscillators. 

Our major aim is to determine the dependence of the transition times between meta- 
stable states, as well as the critical configurations, on the coupling strength, on the whole 
range from weak to strong coupling. This analysis requires a precise knowledge of the sys- 
tem's "potential landscape" , in particular the number and location of its local minima and 
saddles of index 1 [FW981 |Sug96t IKolOni IBE(;tKn4[ IBCiKOSj . In order to obtain this infor- 
mation, we will exploit the symmetry properties of the system, using similar techniques 
as the ones developed in the context of phase oscillators in [AS92t IDGS96a| IDGS96b] . 
for instance. Our study also involves a centre-manifold analysis of the desynchronisation 
bifurcation, which goes beyond existing results on similar bifurcations because a precise 
control of the stability of the bifurcating stationary points is required. 

This paper is organised as follows. Section [2] contains the precise description of the 
model and the statement of all results. After introducing the model and describing its be- 
haviour for weak and strong coupling, in Section [23] we examine the effect of symmetries on 
the bifurcation diagram. A few special cases with small particle number N are illustrated 
in Section [2]6j The desynchronisation bifurcation for general N is discussed in Section [2. 71 
while Section [2.81 considers further bifurcations of the origin. Finally, Section [2.91 presents 
the consequences of these results for the stochastic dynamics of the system. 

The subsequent sections contain the proofs of our results. The proof of synchronisation 
at strong coupling is presented in Section [3l while Section [J] introduces Fourier variables, 
which are used to prove the results for N = 2 and = 3, and for the centre-manifold 
analysis of the desynchronisation bifurcation for general A^. Appendix |A] gives a brief 
description of the analysis of the weak-coupling regime, which uses standard techniques 
from symbolic dynamics, and Appendix [B] contains a short description of the analysis of 
the case = 4. 

The follow-up work [EFGnObj analyses in more detail the behaviour for large particle 
number A^. In that regime, we are able to control the number of stationary points in a 



2 



much larger domain of coupling intensities, including values far from the synchronisation 
threshold. 
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2 Model and Results 
2.1 Definition of the Model 

In our study of the influence of noise on lattice dynamical systems with continuous on-site 
variables, we shall focus on a model which can serve as a paradigm for such systems. 
It is based on coupled bistable ordinary differential equations governed by the compe- 
tition between individual effects and spatial interactions. To name a few examples, it 
applies for instance to chains of particles placed in a quartic potential and coupled by 
springs |FM96j . to stimulus conduction in the myocardial tissue |Kee87] and to certain 
chemical reactions [E N93] . Another motivation for choosing this specific model is that 
it allows for a quantitative description of phenomena, while being generic. In particular, 
it quantifies changes from intensive to extensive properties when the interaction strength 
increases. 

The system is defined by the following ingredients: 

• The periodic one-dimensional lattice is given by A = Z /NX, , where ^ 2 is the 
number of particles. 

• To each site i € A, we attach a real variable Xj G M , describing the position of the ith. 
particle. The configuration space is thus A" = M^. 

• Each particle feels a local bistable potential, given by 

UiO = \c^-lf, eeM. (2.1) 

The local dynamics thus tends to push the particle towards one of the two stable posi- 
tions ^ = lor,^ = — 1. The reason for this choice is that is the simplest possible 
double-well potential, which will be responsible for metastability. Nevertheless, we 
expect this potential to yield a behaviour which is, to a certain extent, generic among 
models in its symmetry class. 

• Neighbouring particles in A are coupled via a discretised-Laplacian interaction, of in- 
tensity 7/2. Such a nearest-neighbour coupling can be expected to yield very different 
dynamics than mean-field models, for instance. 

• Each site is coupled to an independent source of noise, of intensity a. The sources of 
noise are described by independent Brownian motions {Bi{t)}t^Q. 
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The system is thus described by the following set of coupled stochastic differential 
equations, defining a diffusion on X: 

dx^t) = /« (t)) dt + I - 2xf (t) + xUit)] dt + adB,{t) , ieA, (2.2) 

where the local nonlinear drift is given by 

f{C) = -VUiO=C-f ■ (2.3) 
For cj = 0, the system (I2.2p is a gradient system of the form x = —'W^{x), with potential 

V,{x) = U{x,) + ^ j;(x,+i - X,)' . (2.4) 

iGA jGA 

Note that the local potential ?7(^) is invariant under the transformation ^ i— > — ^, im- 
plying that the local dynamics has no preference for either positive or negative ^. An 
interesting question is how the results are affected by adding a symmetry-breaking term 
to C/(^). This question will be the subject of future research. Some preliminary studies 
indicate that several results, such as the presence of synchronisation for strong coupling, 
the structure of the desynchronisation bifurcation, and the qualitative behaviour for weak 
coupling, are not much affected by the asymmetry. However, many details of the bifurca- 
tion diagrams, as well as the metastable timescales and transition paths, will of course be 
quite different. 

2.2 Potential Landscape and Metastability 

The dynamics of the stochastic system depends essentially on the "potential landscape" 
Vy. More precisely, let 

S = 5(7) = {xGX: VV^{x) = 0} (2.5) 

denote the set of stationary points of the potential. A point x G 5 is said to be of 
type (n_,no,n+) if the Hessian matrix of at x has n_ negative, n+ positive and 
riQ = N — n_ — n+ vanishing eigenvalues (counting multiplicity). For each /c = 0, . . . , A^, 
let Sk = <Sk{'y) denote the set of stationary points x G S which are of type {k,0,N — k). 
For k ^ 1, these points are called saddles of index k, or simply k-saddles, while So is the 
set of strict local minima of V^. 

The stochastic system (12. 2p admits an invariant probability measure with density pro- 
portional to e~'^^i(^y'^ J implying that asymptotically, the system spends most of the 
time near the deepest minima of V^. However, the invariant measure does not contain 
any information on the dynamics between these minima, nor on the way the equilibrium 
distribution is approached. Loosely speaking, for small noise intensity a the stochastic 
system behaves in the following way |FW98j : 

• A sample path {x"'{t)}t, starting in a point xq belonging to the deterministic basin of 
attraction A{x*) of a stationary point x* G Sq, will first reach a small neighbourhood 
of X*, in a time close to the time it would take a deterministic solution to do so. 

• During an exponentially long time span, x'^{t) remains in A{x*), spending most of 
that time near x*, but making occasional excursions away from the stationary point. 

• Sooner or later, x"'(t) makes a transition to (the neighbourhood of) another stationary 
point y* £ Sq. During this transition, the sample path is likely to pass close to a saddle 
s € Si, whose unstable manifolds converge to x* and y*. In fact, the whole sample 
path during the transition is likely to remain close to these unstable manifolds. 
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• After a successful transition, the sample path again spends an exponentially long time 
span in the basin of y*, until a similar transition brings it to another point of So (which 
may or may not be different from x*). 

If we ignore the excursions of the sample paths inside domains of attraction, and consider 
only the transitions between local minima of the potential, the stochastic process resembles 
a Markovian jump process on 5o, with exponentially small transition probabilities, the only 
relevant transitions being those between potential minima connected by a 1-saddle. 

Understanding the dynamics for small noise thus essentially requires knowing the graph 
G = {So,£), in which two vertices x*, y* G Sq are connected by an edge e G if and only 
if there is a 1-saddle s G 5i whose unstable manifolds converge to x* and y* . The mean 
transition time from x* to y* is of order e'^^/'^ , where H is the potential difference between 
X* and the lowest saddle leading from x* to y* . 

In our case, the potential V^(.x) being a polynomial of degree 4 in A'^ variables, the set 
S of stationary points admits at most 3^ elements. On the other hand, it is easy to see 
that S always contains at least the three points 

O = (0,...,0), /± = ±(1,...,1) . (2.6) 

Depending on the value of 7, the origin O can be an iV-saddle, or a A;-saddle for any odd 
k. The points always belong to (So, in fact we have for all 7 > 

V^{x)>V^{I+) = V^{r) = -^ Va;G^\ {/-,/+}. (2.7) 

Thus and are the most stable configurations of the system, and also the local minima 
between which transitions take the longest time. 

Among the many questions we can ask for the stochastic system, we shall concentrate 
on the following: 

• How long does the system typically take to make a transition from I~ to /+, and how 

docs the transition time depend on coupling strength 7 and noise intensity cr? 

• How do typical paths for such a transition look like? 

2.3 Weak-Coupling Regime 

In the uncoupled case 7 = 0, we simply have 

5(0) = {-1,0,1}^, |5(0)| = 3^. (2.8) 

Furthermore, 5^(0) is the set of stationary points having exactly k coordinates equal to 
(thus |5fe(0)| = {^^)2^-^). In particular, 5o(0) = {-1, 1}^ has cardinality 2^. 

Hence the graph Q is an A-dimensional hypercube: Two vertices x*, y* G 5o(0) ^'^^ 
connected if and only if they differ in exactly one component. Note that Vq{x*) = —N/4: for 
all local minima x*, and Vq{s) = —{N— l)/4 for all 1-saddles s, implying that all nearest- 
neighbour transitions of the uncoupled system take the same time (of order e^^^"' ) on 
average. 

For small positive coupling intensity < 7 <C 1, the implicit-function theorem guar- 
antees that all stationary points depend analytically on 7, without changing their type. 
In addition, the following result is a direct consequence of standard results on Henon-like 
mappings of the plane: 
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Figure 1. Example of an optimal transition path from to for weak coupling. The 
upper half of the figure shows the local minima and 1-saddles visited during the transition, 
displayed vertically. For instance, the second column means that the first 1-saddle visited 
is a perturbation of order 7 of the stationary point (—1, —1, 0, —1, —1, —1, —1, —1) present 
in absence of coupling. The lower half of the figure shows the value of the potential seen 
along the transition path. 



Proposition 2.1. For any N , there exists a critical coupling 'j*{N) such that for all 
^ 7 < ^*{N), the stationary points x*{j) G 5(7) depend continuously on 7, without 
changing their type. The critical coupling 7*(A^) satisfies 

inf 7*(iV) > \ . (2.9) 

The proof is briefly discussed in Appendix [Aj where we also provide slightly better 
lower bounds on m{N^2^*{N). We expect the critical coupling to be quite close to 1/4, 
however. In particular, we will show below that 7*(2) = 1/3, 7*(3) = 0.2701..., and 
7*(4) = 0.2684... 

Since any stationary point x*(7) = (xf (7), . . . , x'^{'~f)) G ^(7) satisfies x*(7) = x*(0) -|- 
0(7), where each x*(0) is a stationary point of the local potential U, one has 

V,ix*{^))=V,ix*{0)) + Oij') 

N 

= Voix^O)) + J ^{xt^M - xt{0)f + 0(7') . (2.10) 

To first order in 7, the potential's increase due to coupling depends on the number of 
interfaces in the unperturbed configuration x*(0) E 5(0) (recall from (j2.8p that the com- 
ponents of x*(0) only take values in { — 1,0, 1}). The dynamics of the stochastic system 
is thus essentially the one of an Ising-spin system with Glauber dynamics. Starting from 
the configuration I~, the system reaches with equal probability any configuration with 
one positive and — 1 negative spins. Then, however, it is less expensive to switch a spin 
neighbouring the positive one than to switch a far-away spin, which would create more 
interfaces. Thus the optimal transition from I~ to consists in the growth of a "droplet 
of -|- in a sea of — " (Figured]). To first order in the coupling intensity, all visited 1-saddles 
except the first and last one have the same potential value —A/4 -|- 1/4 -|- (8/2)7 -|- 0(7^). 
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Figure 2. (a) Structure of the graph Q in the smaU-coupling regime < 7 < 7*(iV). 
Only edges belonging to optimal paths arc shown: In the Ising-model analogy, these paths 
correspond to the flip of neighbouring spins, (b) The graph Q in the synchronisation 
regime 7 > 71 (A^). 



The energy required for the transition is thus 

v^{{-i, -1, . . . , -1,0, 1,1,...,!) + 0(7)) - v^{r) = J + ^7 + 0(7') , (2.11) 



which is independent of the system size. Note that the situation is different for lattices 
of dimension larger than 1, in which the energy increases with the surface of the droplet 
(cf. [dH04j). 

2.4 Synchronisation Regime 

For strong coupling 7, the situation is completely different than for weak coupling. Indeed, 
we have the following result (see Section [3] for the proof): 



Proposition 2.2. Let 



71 = 71 (^) 



1 A^' 



1 - cos(27r/iV) 27r2 



(2.12) 



Then 5(7) = {0,I^,I } if and only if 'j ^ 71. Moreover, the origin O is a 1-saddle if 
and only if ^ > 71, and in that case its unstable manifold is contained in the diagonal 

V = {x £ X: xi = X2 = ■ ■■ = xn} ■ (2.13) 

We say that the deterministic system is synchronised for 7 > 71, in the sense that the 
diagonal is approached as t ^ 00 for any initial state, meaning that all coordinates Xi are 
asymptotically equal. In other words, the system behaves as if all particles coagulate in 
order to form one large particle of mass A^. The graph Q contains only two vertices 
connected by a single edge (Figure [2)3). By extension, in the stochastic case we will say 
that the system is synchronised whenever all coordinates xf remain close to each other 
most of the time with high probability. Then transitions between /~ and occur in a 
small neighbourhood of the diagonal. 

In this case, the energy required for the transition is 

V,{0)-V,{l-) = ^, (2.14) 
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which is extensive in the system size. Transitions now take a time of order e^l'^'^^ , and are 
thus much less frequent in the synchronisation regime than in the low-couphng regime. 

It is remarkable that as the coupling 7 grows from to 71, the number of stationary 
points decreases from 3"^ to 3. The main purpose of this work is to elucidate in which 
way this transition occurs, and how it affects the transition paths and times. 

2.5 Symmetry Groups 

The deterministic system x = —W^{x) is equivariant (that is, W^{gx) = gWry{x) for 
all X & X) under three different types of symmetries g: 

• Cyclic permutations (corresponding to rotations around the diagonal), generated by 

R{xi, xn) = {X2, XN, xi) , (2.15) 

as a consequence of the particles being identical. 

• Reflection symmetries R'^ S , k = 0, . . . , N — 1, where 

S{XI, . . .,Xn) = {XN,XN-1, . . . ,xi) , (2-16) 

as a consequence of the interaction being isotropic. 

• The inversion 

C(xi, . . . , xn) = -{xi,..., xn) , (2.17) 

as a consequence of the local potential being even. 

The symmetries R and S generate the dihedral group D^, which has order 2N, for N ^ 3, 
and the group Z 2 for iV = 2. For iV ^ 3, the symmetries R, S and C generate a group of 
order 4N, which we shall denote G = Gn = Dn x Z2- For N = 2 the symmetry group is 
the Klein four-group G2 = ^2 x ^2, which as order 4. 

The set of stationary points 5(7), as well as each set Skij) of /c-saddles, are invariant 
under G. Thus G acts as a group of transformations on X, on 5(7), and on each 5/j(7). 
We will use a few concepts from elementary group theory: 

• For X e X, the orbit of x is the set Ox = {gx: g e G}. 

• For X G X, the isotropy group or stabiliser of x is the set Cx = {g & G : gx = x} . 

• The fixed-point set of a subgroup of G is the set Fix(ff) = {x & X: hx = x^h e H}. 
The following facts are well known: 

• For any x, the isotropy group Cx is a subgroup of G and |Ca;||Oa;| = |G| (where |X| 
denotes the cardinality of a finite set X) . 

• For any g E G and x e X, we have Cgx = gCxg~^, so that the isotropy groups of all 

the points of a given orbit are conjugate. 

• For any subgroup H oi G and any g £ G, F[yi(gHg^^) = gF[x{H). 

These facts allow us to limit the study to one point of each orbit, to one subgroup 
in each conjugacy class, and to one type of conjugated fixed-point set. For small N, this 
reduction often suffices to completely determine all stationary points of the system, while 
for larger N, it at least helps to classify the stationary points. 

2.6 Small Lattices 

We now consider some particular cases for illustration. The following applies to the three 
stationary points that are always present: 
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[x2] 
[xl] 

[x2] 
[x4] 



(1,1) 
(0,0) 

(1,-1) 
(1,0) 



1/3 
— I — 



1/2 
— I— 



(x, -x) 



/A 



(0,0) 



7± 
O 





1+ 
O 

r 

FiGURE 3. Bifurcation diagram for the case N = 2 and associated graphs Q. Only one 
stationary point is shown for each orbit of the symmetry group G. The cardinaUty of the 
orbit is shown in square brackets. FuU hues represent focal minima of the potential, while 
dash-dotted lines with k dots represent fc-saddles. 



z* 


0.* 


a* 


Fix(C,0 


(0,0) 


{(0,0)} 


Ga = Z 2 X Z 2 


{(0,0)} 


(1,1) 


{(1,1), (-1,-1)} 


Z2 = {id, 5} 


{(x,x)}:rgK = V 


(1,-1) 


{(1,-1), (-1,1)} 


{id, CS} 


{{x, -x)}^^u 


(1,0) 


{±(1,0),±(0,1)} 


{id} 


{{x,y)}x,ym = X 



Table 1. Stationary points z* G S{0), their orbits, their isotropy groups and the corre- 
sponding fixed-point sets in the case N — 2. 



• For the origin O, Oq = {O}, Co = G and Fix(Co) = {O}. 

• For the global minima I^, Oj+ = Oj- = {/",/+}, Cj± = Dn and Fix(C/±) = V, 
where D is the diagonal defined in (12.131) . 

Case N = 2 

For = 2, we have R = S and the symmetry group is G2 = {id, S,C,CS}. In the 
uncoupled case, the set of stationary points 5(0) is partitioned into four orbits, as shown 
in Table [TJ Figure [3] indicates how the stationary points evolve as the coupling increases 
(the proof is given in Proposition 14. 2p . We see that stationary points keep the same type 
of symmetry as the coupling strength 7 increases, and sometimes merge with a stationary 
point of higher symmetry. 

Below the bifurcation diagram, we show how the corresponding graphs Q change as 
the coupling intensity 7 varies. The two-dimensional hypercube (i.e., the square), present 
for weak coupling, transforms into a graph with two vertices, connected by two edges, 
as the 1-saddles labelled Aa undergo a pitchfork bifurcation at 7 = 1/3 = 7* (2). For 
1/3 < 7 < 1/2, the points with (x, — a;)-symmetry, labelled A, are 1-saddles, representing 
the points of maximal potential height on the two optimal transition paths from I~ to I"*". 
At 7 = 1/2 = 71(2), the 1-saddles undergo another pitchfork bifurcation, this time with 
the origin O, which becomes the only transition gate in the strong-coupling regime. 
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2/3 
— I — 



(2;, X, x) 



(x^x,y)^ 
{x,-x,0) 



B 



■■-9- 

A 
.y i 
/ 

/A 



(0,0,0) 



o 



ix,x,y) 



da 



[x2] (1,1,1) 
[xl] (0,0,0) 
[x6] (0,0,1) 
[x6] (1,-1,0) 
[x6] (1,1,-1) 
[x6] (1,1,0) 

/+ 
O 

I- 

FiGURE 4. Bifurcation diagram for the case iV = 3 and associated graphs G- Only one 
stationary point is shown for each orbit of the symmetry group G. The saddle-node 
bifurcation value is 7* = 7*(3) = (^3 + 2^3 - V3)/3 = 0.2701 . . . 







0.* 


a* 


Fix(a*) 


(0,0,0) 


{(0,0,0)} 


G3 


{(0,0,0)} 


(1,1,1) 


{(1,1,1), (-1,-1,-1)} 


D3 


{{x,X,x)}j:(zM. =V 


(1,-1,0) 


{±(1,-1,0), ±(-1,0,1), ±(0,1,-1)} 


{id, CRS} 


{(x, -a;,0)}:rgK 


(0,0,1) 


{±(0,0,1), ±(0,1,0), ±(1,0,0)} 


{id, RS} 


{{x,X,y)}x,y(^R 


(1,1,-1) 


{±(1,1,-1), ±(1,-1,1), ±(-1,1,1)} 


{id, RS} 


{{x,x,y)}x,ym. 


(1,1,0) 


{±(1,1,0), ±(1,0,1), ±(0,1,1)} 


{id, RS} 


{(a;,a;,y)}a;,ygR 



Table 2. Stationary points z* E 5(0), their orbits, their isotropy groups and the corre- 
sponding fixed-point sets in the case = 3. 



The value of the potential on the bifurcating branches is found to be 



V,{A) 
Vy{Aa) 



>-27)^ ^ 

-^(l-27)'±^(l-37)^ 



(2.18) 



Thus in the case 7 < 1/3 = 7*(2), the typical time for a transition from / to will be 
of order ea+^'yW)/^^'' ^ ^^^^i^ for 1/3 < 7 < 1/2, it is of order e^T^^-^')/'^" , and for 7 > 1/2 
it is of order e^/"^ . 



Case iV = 3 

For = 3, at zero coupling the set 5(0) of stationary points is partitioned into six orbits, 
as shown in Table [2j Their evolution as the coupling increases is shown in Figure U] (for 
a proof, see Proposition I4.3p . The new feature in this case is that two orbits disappear 
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[x2] 



[XI] (0,0,0,0) ^ . 4^i'l'}o 



[x2] (1,-1,1,-1) 
[x4] 



[x4] 



(1,0,-1,0) 
[x4] (1,1,-1,-1) 
[x8] (1,1,0,0) 
[xl6] 
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{x,x,x,x) 



7± 



I 



(x, y.x.y) , , , 
1,0,1,0) ,}-■■---■'■-• — ■■</ / i 
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// 
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{x_,0,_-x_,0) / 

" , y'A 

{x,x,-x,-x^ y 



(x,y,z,t) 



,.yAa 




1+ 
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FlGURE 5. Bifurcation diagram for the case = 4 and associated graphs G (only edges 
corresponding to optimal transition paths are shown). Only one stationary point is shown 
for each orbit of the symmetry group G. The saddle-node bifurcation values are 7* = 
7*(4) = 0.2684 ... and 72 = 0.4004 . . . , while 71 = (3%/2 - 2)/7 = 0.3203 . . . 



in a saddle-node bifurcation at 7 = 7*(3) = (\/ 3 + 2^/3 — \/3)/3 = 0.2701 . . . , instead of 
merging with stationary points of higher symmetry. This accounts for the rather drastic 
transformation of the graph Q from a 3-cube for 7 < 7* (3) to a graph with two vertices 
joined by 6 edges for 7*(3) < 7 < 2/3 = 71 (3). The potential on the ^-branches has value 



2 V 2^' 



(2.19) 



Thus in the case 7 < 2/3, the typical time for a transition from / to will be of order 
g(2+i27-972)/4a2^ while for 7 > 2/3 it is of order e3/2'^'. 



Case iV = 4 

Figure [5] shows the results of a similar analysis for = 4. The determination of this 
bifurcation diagram, which relies partly on the numerical study of the roots of some 
polynomials, is outlined in Appendix [Bj As in the previous cases, a certain number of 
stationary points emerge from the origin as the coupling intensity decreases below the 
synchronisation threshold. In the present case, the desynchronisation bifurcation occurs 
at 7 = 71 (4) = 1, and there are four 1-saddles, labelled A, and four 2-saddles, labelled 
B, emerging from the origin. Two more symmetry-breaking bifurcations affect the A- 
branches, finally resulting in 1-saddles without any symmetry. A second branch, labelled 
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A^'^\ bifurcates from the origin at 7 = 1/2. We do not show the corresponding stationary 
points in the graphs, because they appear always to correspond to non-optimal transition 
paths. 

The examples discussed here give some flavour of how the transition from weak coupling 
to synchronisation occurs in the general case. In the sequel, we will mainly describe 
the desynchronisation bifurcation occurring at 7 = 71 (A^), which can be analysed for 
arbitrary A^. 



2.7 Desynchronisation Bifurcation 

In this section, we consider the behaviour for general particle number N ^ 3, and coupling 
intensities 7 slightly below the synchronisation threshold 71 (A^). Our aim is to obtain 
precise information on 

• the exact number and stability type of stationary points; 

• the potential difference between global minima and saddles of index 1, which governs 
the transition time between these minima; 

• the spatial shape of the critical configurations (i.e., the saddles of index 1), which are 
the highest energy configurations reached during a typical transition. 



Even particle number 

We start by considering the case of even particle number N. The main result is the 
following theorem, which will be proved in Section 14.31 



Theorem 2.3 (Desynchronisation bifurcation, even particle number). Assume that N is 
even. Then there exist S = S{N) > and points A = A{^) and B = -6(7) in X such that 
for 71 — (5 < 7 < 71, the set of stationary points 5(7) has cardinality 2N + 3, and can be 
decomposed as follows: 



So 


= 0j+ 








Si 


= 0a 


= {A,RA,.. 






S2 


= 0b 


= {B,RB,.. 






S3 


= Oo 


= {0}. 







(2.20) 

Furthermore, the value V-f{A) of the potential is the same on all 1-saddles, and is deter- 
mined by 

'-^(1-7)2 ifN = 4:, 

^=1 , ^'-''^ 

[-g (1 - 7/71) +0((1- 7/71)') ifN^6, 

This result tells us that as the coupling intensity 7 decreases below the synchronisation 
threshold 71, there are exactly saddles of index 1 and A^ saddles of index 2 bifurcating 
from the origin. All saddles of the same index belong to the same orbit of the symmetry 
group G. In fact, if the location of one saddle is known, the locations of all other saddles 
of the same type are obtained by cyclic permutation of its coordinates. 

Relation (j2.2ip shows that the height of the 1-saddles decreases as the coupling inten- 
sity 7 decreases. A simple topological argument shows that the unstable manifolds of each 
1-saddle converge to /"*" and /~ . Thus the transition time from I~ to I"*" will be governed 
by the potential difference V^(^) — Vy{I~). 
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-Xi, . . . 


-XL, . . . , -Xi,0) 




4L + 2 


A 


D2 


(xi,... 


; . . . , Xl, - 


-Xl, . . . 


, -Xl+i,...,-Xi) 






B 


D'2 


(xi,... 


, X]^ , X]^ • • • 5 ^1 5 


0, -Xl, 


• • ) ^Lj Xjj, . . . 


-xi,0) 


2L + 1 


A 


(CRS) 


(xi,... 


! ^L; • • • 7 " 


-xi,0) 








B 


{RS) 


(xi, . . . 


, X/^ , Xi^ X\ 


,Xq) 







Table 3. Symmetries of the stationary points bifurcating from the origin at 7 = 71. 
Here L denotes an integer such that N = 4L, = 4L + 2, or TV = 2L -)- 1 depending 
on the value of N (mod 4). The isotropy groups for even N are D2 — {CS, R^^^S) and 
D2 = (CRS, R^/^^^S), where {gi, . . . , g^) denotes the group generated by {51, . . . , gm}- 



N = 4L N = 4L + 2 N = 2L + I 

XL XL ^L+l 




Figure 6. Symmetries of the stationary points A and B bifurcating from the origin at 
7 = 71. Full lines represent symmetry axes, broken lines represent symmetry axes with 
sign change. 



Remark 2.4. The proof of Theorem 12.31 actually yields more precise information on the 
coordinates of the saddles A and B, which can be summarised as follows. 

• The saddles A and B satisfy special symmetries, as shown in Table [3] and Figure [6l 
Namely, their coordinates admit a mirror symmetry, as well as a mirror symmetry 
with sign change, the two symmetry axes being perpendicular. The details depend on 
whether N S 4N or S 4N + 2, because this affects the number of coordinates which 
may lie on the symmetry axes. 

• The value Vy{B) of the potential on the 2-saddles satisfies 

< ^^^""^-^^^^^ ^ - 7/71)^/^) , (2.22) 

indicating that the potential between the two types of saddles becomes flatter as the 
particle number increases. 
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AT G 4N 



W 6 4N + 2 



AT 6 2N + 1 




Figure 7. Components of the stationary points A and B bifurcating from the origin at 
7 = 71, shown for the three different cases N = 4L, N = AL + 2 and N = 2L + 1. 



The components of A = A^j) and B = -6(7) are given by 



Ajij) = Vl-7/71 sm(^ (i - i)) + 0(1 - 7/71) 



5,(7) = -| Vl-7/71 sm(^j) + 0(1 - 7/71) , (2.23) 
(see Figure [7]), except in the case = 4, where 



A{j) = ^/T^ {1,1, -1,-1) and S(7) = yr^(l,0,-1,0) . (2.24) 

For any 7 G [0, 71), there exist stationary points ^(7) and -8(7), satisfying the sym- 
metries indicated in Table [3j The first N/2 components of ^(7) and the first N/2 — 1 
components of B('j) are always positive, and 

lim^(7) = (1, 1, . . . , 1, 1, -1, -1, . . . , -1, -1) , 

lim 5(7) = (1, 1, . . . , 1, 0, -1, -1, . . . , -1, 0) . (2.25) 

We do not claim that ^(7) and -6(7) are continuous in 7 everywhere, though we expect 
them to be so. What we do prove is that for any 7, there is at least one stationary 
point with the appropriate symmetry (and coordinates of the indicated signs). We 
also know that A and B depend continuously on 7 for 7 near and near 71. We 
cannot exclude, however, the presence of saddle-node bifurcations in between. 
We know that the points ^(7) are local minima near 7 = 0. Thus these stationary 
points must undergo at least one secondary bifurcation as 7 decreases. For symmetry 
reasons, we expect that, as in the case = 4 (see Figure [5]), there are two successive 
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symmetry-breaking bifmxations affecting the 1-saddles: First, the mirror symmetry 
with sign change is broken, then the remaining mirror symmetry is broken as well. 
• The error terms in (j2.2ip and (|2.23p may depend on A^. The technique we employ 
here does not allow for an optimal control of the A^-dependence of these error terms 
and of d{N). However, in |BFG06b] . we obtain such a control in the limit of large A^, 
using different techniques. 



Odd particle number 

In the case of odd particle number N, the results are slightly weaker, since we are not able 
to obtain such a precise control on the number of 1-saddles and 2-saddles created in the 
desynchronisation bifurcation. The main result is the following theorem, which will also 
be proved in Section 14. 3i 



Theorem 2.5 (Desynchronisation bifurcation, odd particle number). Assume that N is 
odd. Then there exists 6 = S{N) > such that for 71 — 5 < 7 < 71, the set of stationary 
points 5(7) has cardinality MN + 2>, for some I 1. All stationary points are saddles of 
index 0,1,2 or 3, with 

So = 0j+ ={!+,!-} , 

Ss = Oo = {O} , (2.26) 

and \Si\ = \S2\ = 2iN . The value V-y{A) of the potential on any saddle A of index 1 
satisfies 

^ = -^(1 - 7/71)' + 0{{1- 7/71)') . (2.27) 

In fact, we expect that ^ = 1, so that there are exactly 2A^ saddles of index 1 and 2A^ 
saddles of index 2. This would lead to the following conjecture. 

Conjecture 2.6 (Desynchronisation bifurcation, odd particle number). For odd N and 

7i — (5 < 7 < 71, the set of stationary points 5(7) has cardinality 4:N + 3. There exist 
points A = A{^) and B = -8(7) in X such that S can he decomposed as follows: 



So 


= Oj+ 








Si 


= Oa 


= {A,RA,... 


,R^-^A,-A,-RA,... 




S2 


= Ob 


= {B,RB,.. 


. ,R^~^B,-B,-RB,.. 


.,-R''-'B}, 


S3 


= Oo 


= {0}. 







(2.28) 

This conjecture would follow as a consequence of a much simpler conjecture on the 
behaviour of certain coefficients in a centre-manifold expansion, which can be computed 
iteratively, see Section 14. 3i We know that it is true for = 3, and it can be checked by 
direct computation for the first few values of A^. Numerically, we checked the validity of 
the conjecture for A^ up to 101. In |BFG06b] . we show that the conjecture is also true for 
sufficiently large A^. 

Remark 2.7. Again, the proof of the theorem also yields more precise information on the 
location of the saddles. 



• The potential difference between the 2-saddles and 1-saddles satisfies 

0^}m_}m^O{il-jM-). (2.29) 
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• If the conjecture is true, then the components of ^ = ^(7) and B = B(j) are of the 
form 



^i(7) = ^ Vl-7/71 ^'""[^j) + 0(1 - 7/71) , 

5,(7) = ^ - 7/71 cos(^j) + 0(1 - 7/71) , (2.30) 

and satisfy the symmetries indicated in Table [3l 

• Independently of the validity of the conjecture, for any 7 G [0, 71), there exist station- 
ary points ^(7) and -8(7), satisfying the symmetries indicated in Table [3l The first 
(A^ — l)/2 components of ^(7) are strictly positive, and 

lim^(7) = (1,..., 1,-1,..., -1,0) . (2.31) 
7— »o 

2.8 Subsequent Bifurcations 

We shall show in Section [3] that the origin undergoes further bifurcations at 

in which the index of the origin O increases by 2 (except for the case where N is even 
and M = N/2, where the index increases by 1), and new saddles A^^'^ and B^^'^ of 
index 2M — 1 and 2M are created. Consequently these saddles are not important for 
the stochastic dynamics, and we shall not provide a detailed analysis here. We briefly 
mention a few properties of these bifurcations, which we will prove in |BFG06b] to hold 
for sufficiently large N/M: 

• The number of newly created stationary points is given by 4A^/ gcd(A'^, 2M), where 
gcd(iV, 2M) denotes the greatest common divisor of N and 2M. 

• The number of sign changes of function of j for these new stationary points is 
equal to 2M. M can therefore be considered as a winding number. 

• If iV and M are coprime, the new stationary points x satisfy the symmetries shown 
in Table El while for other M they belong to larger isotropy subgroups. 

Example 2.8. For = 8, the origin bifurcates four times as 7 decreases from -|-cxd to 
0. We show the symmetries of the bifurcating stationary points in Table HI They are 
obtained in the following way: 

• Compute the eigenvectors of the Hessian of the potential at the origin; 

• Determine the corresponding isotropy subgroups of Gg; 

• Write the equation i = —W^{z) restricted to the fixed-point set of each isotropy 
subgroup, and study the bifurcations of the origin in each restricted system. 

For instance, for winding number M = 1 and orbits of type A, we obtain 

X = (1- |7)x + ^-fu - x^ , 

1 s (2-33) 

y = ^'yx + 27)y - y ■ 

The origin bifurcates for 7 = 2 it \/2- An analysis of the linearised system shows that 
for 7 slightly smaller than 2 -|- \/2, the new stationary points must lie in the quadrants 
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3 
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-x,0) 


4 
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8 


(x, 


-X, 






—x) 







Table 4. Fixed-point sets of the stationary points bifurcating from the origin for = 8, 
for different winding numbers M. Points of winding number M — I and M — 3 actuaUy 
have the same fixed point spaces, but we choose the signs of the components in such a way 
that X and y always have the same sign for the actual stationary points. 



{(x,y): X > 0,y > 0} and {{x,y): x < 0,y < 0}. These quadrants, however, are invariant 
under the flow of x = — VV^(x), since, e.g., x > if x = and y > 0, and y > if y = 
and X > 0. Hence the points created in the bifurcation remain in these quadrants. The 
points created at 7 = 2 — ^/2, which correspond to the winding number M = 3, he in the 
complementary quadrants {(x, y) : x > 0, y < 0} and {(x, y) : x > 0, y < 0}. 
For 7 = 0, the only points with the appropriate symmetry are 

A«(0) = (1,1,1,1,-1,-1,-1,-1) , S«(0) = (1,1,1,0,-1,-1,-1,0) , 
^(2)(o) = (1,1,-1,-1,1,1,-1,-1) , S(2)(0) = (1,0,-1,0,1,0,-1,0) , 

^(3)(0) = (1,-1,-1,1,-1,1,1,-1) , ^(i^_;l^i^o^_;^^1^_1^q) ^ 

^(^)(0) = (1,-1,1,-1,1,-1,1,-1) . (2.34) 

Note that the cases M = 2 and M = 4 are obtained by concatenation of multiple copies 
of stationary points existing for = 4 and N = 2, respectively. 

2.9 Stochastic Case 

We return now to the behaviour of the system of stochastic differential equations 

dx^(t) =/(xf(t))dt + |[xf+i(t)-2xf(t) + x^_i(t)]dt + ad5i(t) , ieA. (2.35) 

Recall that our main goal is to characterise the noise-induced transition from the con- 
figuration /~ = (—1, —1, . . . , —1) to the configuration /+ = (1, 1, . . . , 1). In particular, 
we are interested in the time needed for this transition to occur, and in the shape of the 
critical configuration, i.e., the configuration of highest energy reached in the course of the 
transition. 

Since the probability of a stochastic process in continuous space hitting a given point is 
typically zero, we have to work with small neighbourhoods of the relevant configurations. 
Given a Borel set A C X, and an initial condition xq ^ X \ A, we denote by t^^^{A) the 
first-hitting time of A 

T^'\A) = mf{t > 0: x^{t) G A} . (2.36) 
Similarly, for an initial condition xq € A, we denote by t^"^{A) the first-exit time from A 

^cxit^^) = inf{t > 0: x'^(t) ^ ^} . (2.37) 

We can now formulate our main results, which are similar in spirit to Theorems 3.2.1 
and 4.2.1 of [dH04j . 
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2/3 1 7/7i(2) 1 7/7i(3) 



Figure 8. The normalised potential differences h^ij) for = 2 and = 3. Tlie broken 
curve for iV = 2 sliows tlie potential difference {Vj{A) — V^(/^))/2, wfiich is smaller than 
hN{l) for 7/71 < 2/3 (compare (|2.18l) ). because for these parameter values, the 1-saddle 
is the point labelled Aa. For general particle number iV, we know that /ijv(7) behaves like 
1/4 - cn{1 - ihif as 7 / 71, and that /lAr(O) = 1/47V. 



Theorem 2.9 (Stochastic case, synchronisation regime). Assume that the coupling strength 
satisfies 7 > 71 = (1 — cos(27r/A^))~"'^ . We fix radii < r < R < 1/2, and denote by 
= T^^^{B{I^,r)) the first-hitting time of a hall B{I^,r) of radius r around . Then 
for any initial condition xq € B{I^,r), any N ^ 2 and any 6 > 0, 



and 



lim P^o {eW^-^y^" < r+ < eW^+^)/-' } = 1 (2.38) 



lim log E'^o {r+} = — . (2.39) 



Furthermore, let tq = T^^^{B{0,r)), and let 

r_ = inf{t > r^^'*(i3(/-,i2)): xt G B{I- ,r)} (2.40) 

he the time of first return to the small hall B{I~ ,r) around after leaving the larger hall 
B{I-,R). Then 

lim P^ojro < r+ I T+ < r_} = 1 . (2.41) 



We will omit the proof of this result, which is a direct consequence of Proposition 
and standard results in Wentzell-Preidlin theory, see for instance |FW98t |Sug96[ IKif81| . 

Relations ()2.38p and (j2.39p mean that the transition from /~ to typically takes a 
time of order e^/'^'^ . Relation ()2.4ip shows that, provided a transition from B{I^,r) to 
B{I^,r) is observed, the process is likely to pass close to the saddle at the origin on its 
way from one potential well to the other one. The origin thus plays the role of the critical 
configuration in the synchronisation regime. 

Theorem 2.10 (Stochastic case, desynchronised regime). Assume 7 < 71. Let r,R and 
and r_ be defined as in Theorem \2.9[ and fix an initial condition xq G B{I~,r). Let 
hnil) = (y-yi^) ~^'y{L^)) denote the normalised potential difference between 1-saddles 
and potential minima at L^ . Then /iAr(7) satisfies 

hN{l) = \-CN{l-lhif + 0{{l-^hif) as^/^i, (2.42) 
where C2 = C4 = 1/4 and cn = 1/6 for N = 3 and all N ^ 5, and one has 

lim^P"°{e(2^''^W-^)/'^' < r+ < e^^Nhr,{^)+S)/a^ ^ ^ ^ ^ ^2.43) 
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and 

lim a"^ logE^" {r+} = 2NhN{'~i) . (2.44) 
Furthermore, assume that either N is even, or N is odd and Conjecture \2.6\ holds. Let 

TA = T'^'{\jB{gA{^),r)) . (2.45) 

geG 

Then there exists a 5 = d{N) > such that for 71 — 5 < 7 < 71 

lim P^ojrA < T+ I r+ < T_| = 1 . (2.46) 

The proof is again standard, and we will omit it here. 

Relations (|2.43p and (j2.44p mean that the transition from /~ to J"*" typically takes a 
time of order e^^^Nii)/'^ ^ while relation (j2.46p shows that the set of critical configurations 
is the orbit of ^(7). 

We conclude the statement of results with a few comments: 

• The very precise results in [BEGK04t IBGK05) allow in principle for a more precise 
control of the expected transition time r+ than the exponential asymptotics given 
in (|2.39p and (12.441) . However, these results cannot be applied directly to our case, 
because they assume some non-degeneracy condition to hold for the potential, which 
excludes symmetries as in the system we are investigating. Therefore, in the present 
work we content ourselves with the exponential asymptotics. A finer analysis is feasi- 
ble, and we plan to provide such a further study. 

• We have seen that the potential difference between the 1-saddles A and the 2-saddles 
B becomes smaller and smaller as the particle number increases. As a consequence, 
for any given noise intensity o" > 0, the transition paths become less localised as 
the particle number A*" increases. This reflects the fact that the system becomes 
translation- invariant in the large- A^ limit. 

• As the coupling intensity 7 decreases, the configurations of highest energy reached 
in the course of a transition become more and more inhomogeneous along the chain, 
and we expect them to converge to configurations of the form (1, . . . , 1, 0, —1, . . . , —1) 
in the uncoupled limit. In addition, several local minima and saddles should appear 
by saddle-node bifurcations along the optimal transition path, thereby increasing the 
number of metastable states of the system. 
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3 Lyapunov Functions and Synchronisation 

We now turn to the proofs of the statements made in Section [21 and start by introducing 
a few notations. The "interaction part" of the potential is proportional to 

Wix) = ^ Y^ixi - x^+lf = l\\x- Rxf = {x, Ex) , (3.1) 

where S is the symmetric matrix S = 1 — ^(i? + R^^). Hence, the potential Vy{x) can be 
written as 

Vy{x) = (1 - 7S)x) + ^ x^ (3.2) 

The eigenvectors of S are of the form Vk = (1, w'^, . . . , u;^^"^^'^)'^, A; = 0, . . . , — 1, where 
CO = e^'^^l^ , with eigenvalues 1 — cos{2TTk/N). This implies in particular that the Hessian 
of the potential at the origin, which is given by 7S — 1, has eigenvalues — A^, where 

/ 2'Kk\ 7 , , 

Afc = A_fe = 1 - 7 1 - cos — - =1 . (3.3) 

V N J 7fc 

The origin is a 1-saddle for 7 > 71. As 7 decreases, the index of the origin increases by 2 
each time 7 crosses one of the 7^, until it becomes an A^-saddle at 7 = 7LAr/2j • 

We now show that for 7 > 71 (i.e., Ai < 0), W{x) is a Lyapunov function for the 
deterministic system x = — Wy{x). 

Proposition 3.1. For any initial condition xq, the solution x{t) of x = —'Wy{x) satisfies 
^W{xit)) ^ 2(1 - 7/7i)Ty(x(t)) - ^W{xit)f . (3.4) 

As a consequence, if ^ is strictly larger than 71, then x{t) converges exponentially fast to 
the diagonal, and thus the only equilibrium points of the system are O and I^. 

Proof: We first observe that the relation 

/(xi) - /(xi+i) = {xi - Xi+i) [1 - (xf + XiXi+i + xf+i)] (3.5) 

allows us to write 

— (x-Rx) = U{x,-f){x-Rx), (3.6) 
at 

where n(x, 7) = 1 — 7S — D{x). Here D[x) is a diagonal matrix, whose ith. entry is given 
by xf + XiXi+i + x'jj^^, and can be bounded below by \{xi — Xj+i)^. It follows 

^W{x{t)) = {x- Rx, ^{x - Rx)) 

= (x — -Rx, n(x, 7)(x — Rx)) 

^ (x - Rx, (1 - 7S)(x - Rx)) - ^ ^{xi - Xi+i)^ 

ieA 

^ Ai||x-i?xf - -47||x-i?x||^ (3.7) 

ill II ^^11 II V J 

by Cauchy- Schwartz. This implies (|3.4p . If 7 ^ 71, then W{x(t)) converges to zero as 
t — > 00 for all initial conditions, which implies that all stationary points x* must satisfy 
W{x*) = 0, and thus lie on the diagonal. The only stationary points on the diagonal, 
however, are O and /±. □ 
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Proof of Proposition 12.21 The assertion on the unstable manifold follows from the 
invariance of the diagonal and the fact that O is a 1-saddle if and only if 7 > 71. For 
7 < 7i, we will show independently that there exist stationary points outside the diagonal. 
Note however that relation (j3.4p shows that these points must lie in a small neighbourhood 
of the diagonal for 7 sufficiently close to 71. □ 

Let us also point out that the growth of the potential away from the diagonal can be 
controlled in the following way. 

Proposition 3.2. For any G T> and xj_ orthogonal to the diagonal T>, the potential 
satisfies 

V^{xu+x^) >V.,{xu) + l(^-l)\\x^f . (3.8) 
' " 2 V71 / 

Proof: Using p.2p and the fact that T,x\\ = 0, we obtain for any A G M 

+ = -2(11^11 II' + -^'ii^^ii') + 4 Ed^ii + ^^J^)^ + ^^^) ■ (3-9) 

i=l 

The scalar product {x_l,T,x±) can be bounded below by ||x_l|P/7i. Moreover, applying 
Taylor's formula to second order in A, and using the fact that all components of xy are 
equal while the sum of the components of x± vanishes, (j3.8|) follows. □ 



4 Fourier Representation 

Let uj = e^'^ ' . The Fourier variables are defined by the linear transformation 

yk = ^Yl ^'''''1 ' keA* = Z/NZ . (4.1) 

The inverse transformation is given by 

Xj = ^''Vk , (4.2) 
fceA* 

as a consequence of the fact that Ylf=i^''^''~^'^ — ^^ki- Note that = y^, so that we 
might use the real and imaginary parts of yk as dynamical variables, instead of yk and 
y-k- The following result is obtained by a direct computation. 

Proposition 4.1. In Fourier variables, the equation of motion x = — VKy(x) takes the 
form 

Vk = ^kVk - ^ ykiVkiVka , (4.3) 

fei,A;2,fc3GA* 
ki+k2+k3=k 

where the A^ are those defined in (|3.3p . Furthermore, the potential is given in terms of 
Fourier variables by 

%{y) = -y X] ^k\yk?' + J X] VkiVk^yksyk^ ■ (4.4) 

fceA* fci,A:2,fe3,fc4GA* 

fcl+fc2+fc3+fe4=0 
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R 


Xj 


^ Xj + 1 


Vk 




RS = SR ^ 


Xj 




Vk 


^yk = y-k 


C 


Xj 


-Xj 


Vk 


^ -Vk 



Table 5 . Effect of a set of generators of the symmetry group Gn on Fourier variables. 

The effect of the symmetries on the Fourier variables is fully determined by the action 
of three generators of the symmetry group G, as shown in Table [5j 

A particular advantage of the Fourier representation is that certain invariant sets of 
phase space take a simple form in these variables. For instance, for any £, 

xe_j = Xj \/j => Lj'^^yj: = ykyk yk = J^''^rk V/c (4.5) 

where the are all real. Similarly, for any ^, 

xii-j = -Xj Mj uj^^Vk = -Vk V/c 

where the are all real. 

4.1 The Case N = 2 

For = 2, we have u; = — 1, and the Fourier variables are simply i/q = (xi + X2)/2, 
yi = {x2 — xi)/2. The equations (|4.3p become 

yo = yo[l- {yl + 3yf )] , 
yi = yi[Ai - (3^0 + y?)] > 

with Ai = 1 — 27. The potential is given by 

%{y) = -yl - Aiyf + \{yt + Qyhl + yi) . (4.8) 

Proposition 4.2. The bifurcation diagram for N = 2 is the one given in Figure\^ and 
the value of the potential on the bifurcating branches is given by 

V,iA) = -hj , V,{Aa) = ^(A? - 6A1 + 1) . (4.9) 

Proof: In addition to the origin, there can be three types of stationary points: 

• If yi = 0, yo / 0, then necessarily 7/0 = il) yielding the stationary points in 
original variables. 

• If yo = 0) yi 7^ 0, there are two additional points A, RA, given by yi = ±-v/Ai whenever 
Ai > 0, i.e., 7 < 1/2. In original variables, these have the expression (it-v/Ai, =FV^)) 
so that they have the {x, — x)-symmetry. 

• If yo, yi 7^ 0, there are four additional points, given by 8yg = 3Ai — 1, Sy\ = 3 — Ai, 
provided Ai > 1/3, i.e., 7 < 1/3. 

It is straightforward to check the stability of these stationary points from the Jacobian 
matrix of (|4.7p . and to compute the value of the potential, using (14. Sp . □ 



^ yfc = ia;^^/VfcVA: (4.6) 
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4.2 The Case = 3 

For = 3, we choose A* = {— 1,0,1}. The equations in Fourier variables read 

yo = yo-(yo + y? + yr^ + 6yo|yiP) , 



yi = Aiyi - 3(yi|yi| + yoyi + yoVi) , 



with Ai = 1 - I7. 



Proposition 4.3. The bifurcation diagram for N = 3 is the one given in Figure^ where 
the saddle-node bifurcations occur for 



7 = 7^(3) = ^3 + 2^-^ ^ Q 2701 .... (4.11) 
3 

On the 1- saddles, the potential has value 

V,{A) = -Ia? . (4.12) 
Proof: Using polar coordinates yi = ri e"^^, the equations (|4.10p become 

yo = yo(i - yo - ^^i) - 2?"i cos 3ipi , 

h = n [Ai - 3{rl + yori cos 3(/?i + yg)] , (4-13) 
^1 = ^yon sin 3931 . 

In addition to the origin, there can be three types of stationary points: 
• If ri = 0, yo 7^ 0, then necessarily j/q = il; yielding the points 



• If yo = 0, ri 7^ 0, we obtain six stationary points given by ri = yAi/S and cos 3(/?i = 0, 
provided Ai > 0, that is, 7 < 2/3. These points have one of the symmetries {x, —x, 0), 
(x,0, — x) or (0, X,— x). 

• If sin3(/9i = 0, it is sufficient by symmetry to consider the case tpi = (i.e., yi real). 
These points have the (x, x, y)-symmetry. Setting yo = u + v, ri = u — v, we find that 
stationary points should satisfy the relations 

Ai/3 = 3u^ + v'^ , 

, . (4.14) 

= 24ii2t; + (1 - AOii + (1 - |Ai)^; . 

Taking the square of the second equation and eliminating u yields a cubic equation 
for v'^. In fact, the variable z = v'^ — {1 + Ai)/12 satisfies 

z-Xz + ii = 0, ^^48' 1728 ' ^ ' 

This equation has three roots for Ai slightly smaller than 1, and one root for Ai = 0. 
Bifurcations occur whenever the condition 27//^ = 4A^ is fulfilled, which turns out to 
be equivalent to (1 — Ai)^g(Ai) = 0, where 

5(Ai) = A\\ - 16A? + 12A? - 4Ai + 1 . (4.16) 

Since 5(0) = 1 and g[l) = —3, and it is easy to check that < on [0, 1], there can be 
only one bifurcation point in this interval, whose explicit value leads to the bifurcation 
value dm]). □ 
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4.3 Centre-Manifold Analysis of the Desynchronisation Bifurcation 

Assume N 3. We consider now the behaviour for 7 close to 71, i.e., for Ai close to 0. 
Setting z = (uq, y2,y-2, ■ ■ ■), the equation (|4.3p in Fourier variables is of the form 

ilk = hyk + 9k{yi,yi,z) , (4.17) 

where 

9k{yi,yi,z) = - ^ yk^ykiyk-j ■ (4.18) 

ki+k2+k3=k 

For small Ai, the system admits an invariant centre manifold of equation 

Vk = hk{yi,mAi) , k^±l, (4.19) 
where the /i^ satisfy the partial differential equations 



>^khk + gk{yi,yiAf^j}j) = ^[-^^^i +5i(yi,yi,{^i}j)] + ^[Aiyi + gi{yi,yi,{hj}j)] . 

(4.20) 

We fix a cut-off order K. For our purposes, K = 2N will be sufficient. We are looking for 
an expansion of the form 

hk{yi,m,Xi)= /iL(Ai)yi"yr'" + 0(|yi|^). (4.2i) 

n,m^O 

First it is useful to examine the effect of symmetries on the coefficients. 

Lemma 4.4. The coefficients in the expansion of the centre manifold satisfy 
. /iJ„(Ai) G R^^ 

• ^nm(Al) = /^mn (Al); 

• ^nm(Ai) = ifn — mj^k (mod A^); 

• /i^^(Ai) =0 if n + m is even. 

Proof: The centre manifold has the same symmetries as the equations (j4.3p . Thus 

• The i2-symmetry requires u:i^hk{yi,yi,\i) = hk{u!yi,uJyi, Xi), yielding the condition 
{lu^ — u;"~™')/iJ;^(Ai) = so that h^^{Xi) vanishes unless n — m = k (mod iV); 

• The iJS-symmetry requires /ifcd/ii yTi Ai) = /ifc(yr, yi, Ai), and yields the reality of the 
coefficients; 

• The C-symmetry requires — /ifc(yi, yl, Ai) = /ifc(— yi, — yl, Ai), yielding the condition 
((—1)"^™ _l_ l)/i^^(Ai) = 0, so that /i^„(Ai) vanishes unless n + m is odd; 

• The condition /i.-fc(yi, yl, Ai) = /ifc(yi) yT) Ai) yields the symmetry under permutation 
of n and m. □ 



From now on, we will write n — m = k instead of n — m = k (mod N). Lemma W. 
allows us to simplify the notation, setting 

hkiyi,m,Xi)= Yl /in™(Ai)yi"yr"' + 0(|yi|'^) . (4.22) 

n—m=k 
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It is convenient to set yT, Ai) = yi and yl, Ai) = yT. Then (j4.22p holds for 

/c = ±1 as well if we set hnm = ^ni^mo whenever n — m = 1, and hnm = SnO^mi whenever 
n — m = —1. The equation on the centre manifold can be written as 

yi = Aiyi - ^ /ifci(yi,yr, Ai)/ifc2(2/i,yr, Ai)/ifc3(yi,yr, Ai) 

fcl+fe2+fc3 = l 

= Aiyi - Cnm{\i)Viyr + 0{\yi\'') , (4.23) 



n— m=l 



where 



Cnm (Al ) — ^2 /inimi(Al)/ln2m2(Al)^n3m3(Al) G M . (4.24) 

rii^O: ni+n2+n3=n 
mi^O: mi+m2+m3=m 

Note that c„m(Ai) = whenever n + m is even. In polar coordinates yi = rie^'^'^, 
Equation (|4.23|) becomes 

ri = Airi- YI c„^(Ai)r]^+"cos((n- m- l)<^i) +0(rf) , 

^i = - Y c„„(Ai)rl^+'"-i sin((n - m - l)<^i) + 0(rf -^) . 

n,m^O, 3^n+m<ii' 
n— m=l 

In general, C21 = S/ifg/ioi = 3 is the only term contributing to the third-order term of 
ri. The only exception is the case = 4, in which C03 = hf^ = 1 also contributes to the 
lowest order, yielding 

[Am - (3 + cos499i)rf + C'(rf) if iV = 4 . 

This shows that all stationary points bifurcating from the origin lie at a distance of order 
\fXi from it: They satisfy ri = y^Ai/(3 + cos 4991) + 0{Xi^'^) in the case = 4, and 

ri = -s/ Ai/3 + ©(A^^^) otherwise. 

The terms with n — m = 1 do not contribute to the angular derivative ipi. In the 
particular case = 4, we have 

ipi = sin(4(pi)r? + 0{rf) , (4.27) 

yielding 8 stationary points, of alternating stability. Otherwise, we have to distinguish 
between two cases: 

• If is even, the lowest-order coefficient contributing to 921 is co,Ar-i, giving 

= co,7V-i(Ai)rf-2 sm{N^i) + 0{r^~^) . (4.28) 

Thus if we prove that co,Af-i(Ai) / 0, we will have obtained the existence of exactly 
2N stationary points, of alternating stability, bifurcating from the origin. 
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• If is odd, then n + m is even whenever n — m = ibA^ + 1, which imphes by Lemma [4.4l 
that Cnm = for these {n,m). The lowest-order coefficient contributing to ipi is thus 
Co,2Ar-i, giving 

'Pi = co,2JV-i(Ai)r2^-' sm{2N^i) + ©(r^^-^) . (4.29) 

Thus if we prove that co,2Af-i(Ai) ^ 0, we wih have obtained the existence of exactly 
4N stationary points, of alternating stability, bifurcating from the origin. 

Remark 4.5. Let ri and r[ be solutions of the equation ri = obtained for two different 
values of (pi, say pi = and ipi = n/N. For even A'^, using the fact that the terms 
up to order — 2 in the first equation in (|4.25p do not depend on ipi, one can see that 
r[ — ri = 0{x\^ ^^^^). For odd A^, one obtains in a similar way r[ — ri = 0{X^i^ ^^^'^). 

In order to compute sign of the coefficients cq^n-i or co,2Af-i) we need at least to know 
the coefficients /lom foi' odd m up to A^ — 3 or 2A^ — 3, respectively. By continuity, however, 
it is sufficient to compute them for Ai = 0. We henceforth set hnm = hnm{0). 

Lemma 4.6. For all odd m > such that m ^ ±1, /iQm = hmo satisfies 

Am^Om — ^ ^ ^0mi^0m2 ^07713 

TUi^O: nii+ni2+m3=m 

v^O: v=m+l 
rrii^O: mi+m2+m3+v='m 



- ^ vhnj^ohn^ohnaohov ■ (4.30) 

v>0: v=m 
rti^O: ni+n2+n3+v=m+l 

Furthermore, if either N is even and 1 ^ m ^ N — 3, or N is odd and 1 ^ m ^ 2A^ — 3, 
then 

Am^Om = ^0mi^0m2^0m3 ■ (4-31) 

rrii^O: m-i+m2+m2=m 

Proof: By invariance of the centre manifold, hk{yi,yT, 0) has to satisfy the equation 

dhk , _ r, . . , dhk- 



Afc/ifc = -gk{yi,yi,{hj}j) + -^9i{yi,yi,{hj}j) + -^giiyi,yi,{hj}j) ■ (4.32) 

oyi oyi 
Plugging in the series (j4.18p of gk and (j4.22p of /i^, this can be seen to be equivalent to 

^n—m^nm — ^ ^ hn-^rn\hn2m2^n^'m^ 

rii^O; ni+n2+n3=n 
TTii^O: mi+m2+m3=m 

^ ^ ^^nimi ^n2m2 ^".3(713 ^nt) 

u,v'^0: u—v=n—m 
rii^O: ni+n2+n3+u=n+l 
rrii^O: mi+m2+m3+v=m 



^ ^ 'y^mmi ^n2r7i2 ^nsms ^nt) • (4.33) 



u,v^O: u—v=n—m 
rii^O: ni+n2+n3+v=m+l 
nii^O: mi+m2+m3+u=n 
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In the special case n = 0, the second sum vanishes unless u = 1, and ni = n2 = = 0. 
In the third sum, we must have v > and u = mi = m2 = = 0. Finally, the fact that 
A-m = A^ yields gJO]). 

Assume now that is even and m ^ N — 2>. In the second sum, v cannot exceed 
m, but then the condition v = m + 1 would require m ^ N — 1. Thus the second sum 
vanishes. In the third sum, v cannot exceed m + 1, and thus v = m. However, in that 
case rai + n2 + n3 = 1, so that two must be zero. Since /iqo = 0, the third sum vanishes 
as well. 

If is odd and m ^ 2N — 3, then the second sum in (j4.30p allows for v = m + 1 — N . 
Then, however, we would have mi + m2 + m^ = N — 1, which is even. Thus at least one of 
the mi is even, yielding a vanishing summand. The third sum in (|4.30p allows for v = m 
and V = m — N . In the first case, however, the summand vanishes for the same reason 
as before, while in the second case, we would have rii + n2 + = N + 1, which is even. 
Thus at least one of the rii is even, yielding again a vanishing summand. □ 

Proposition 4.7. If N is even, then 



As a consequence, for Ai > sufficiently small, the system admits exactly 2N stationary 
points on the centre manifold. The points with ipi = 2k7r/N have one stable and one 
unstable direction if N £ 4N , and two stable directions if N € 4N + 2, and vice versa for 
the points with ipi = {2k + 1)it/N . 

Proof: By (j4.24p . it is sufficient to compute /lom = ^mo ^or odd m between 1 and A^ — 3. 
Recall that Aq = 1, and A^ = A_fc < for 3 ^ A; ^ A^ — 2. Using Hqi = 1 as starting point, 
it is easy to show by induction that sign(/io,2£+i) = because all summands have the 

same sign at each iteration. Likewise, all summands of cq at-i have sign (-1)^/2. □ 

Proof of Theorem 12.31 We consider the case A^ = 4L, the proof being similar for 



• First note that the set B = {y: yt £ ^ VA:} is invariant under the dynamics (it 
corresponds to Vj). The intersection of B with the centre manifold is 

one-dimensional, and can be parametrised by ri G M (while ipi =0). Since ri = 
Airi — 3rf + 0{rf), the system admits at least three stationary points O and ±5' in B 
for small positive Ai. The stationary point B' is stable in the ri-direction, and unstable 
in the (/9i-direction. Since there are A^ — 3 stable and 1 unstable directions transversal 
to the centre manifold, B' is a 2-saddle. The same holds for the cyclic permutations 
RB', . . . ,R^~^B', which correspond to (pk = 2k7r/N, and thus lie in the fixed-point 
sets of conjugate symmetry groups. Applying the inverse Fourier transformation ()4.2p . 
we find that the coordinates of B' in X satisfy 



Setting B = R-^/^B' yields the expression ()2.23p for the coordinates. 
• A similar argument shows the existence of A^ stationary points A,RA,...,R^~^A, 
corresponding to = {2k + l)7r/N, which are 1-saddles because they are stable in 
the (/^i-direction. In X, their coordinates satisfy one of the symmetries Xn^-j = —Xj, 
no = l,...,N. 




(4.34) 



N = 4L + 2. 




(4.35) 
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• For Ai small enough, the equation (fi = admits exactly 2N solutions, so that there 
are no further stationary points on the centre manifold. Proposition 13.11 shows that 
for any ry > 0, we can find a 6 > such that if ^ > ji — 5, there can be no stationary 
points outside an //-neighbourhood of the diagonal. Together with a local analysis near 
the diagonal, and the fact that the centre manifold is locally repulsive, this proves that 
there are exactly 2N + 3 stationary points. 

• Consider now the set 

A+ = {x e X: Xj = -XN+i-j = XN/2+i-j Vj, xi,. . . ,XL > 0} . (4.36) 

We claim that this set is positively invariant under the flow of x = — VV^(x). Without 
the condition xi, . . . ,xl > 0, the invariance follows from equivariance. Now if xj = 
for some j while the other Xi are positive, one easily sees that —dx^V^ = {xj^i + 
Xj+i)7/2 is positive, showing the invariance of ^4.. Since the potential is increasing 
at infinity, there must be at least one stationary point in A-^-, which we denote ^(7). 
As 7 ^ 0, the only stationary point in A+ is the point (1, . . . , 1, —1, . . . , —1). We 
proceed similarly for B['^). 

• Finally, the value of the potential at the stationary points can be computed with the 
help of the expression (j4.4p for the potential in Fourier variables. The only terms 
contributing to leading order in Ai are the term Ai|yip in the first sum, and terms 
of the form yfy'^i in the second sum. The remaining terms are of smaller order. The 
relation on the difference Vy{B) — V^{A) is a consequence of Remark 14.51 This proves 
the theorem. □ 

In the case of odd N, the situation is more difficult, because not all summands in the 
recursion defining the hom are of the same sign. A partial result is 

Lemma 4.8. // N is odd, then 

sign(V2£+i) = (-1)' /or£ = 0,l,...,^^-l . (4.37) 

Proof: As before, using the fact that A^ = A_fc < for 1 ^ ^ iV — 1. □ 

Conjecture 4.9. For £ = {N - l)/2, N - 2, hom has sign (-l)'^+^ and 

co,2N-i > . (4.38) 

Numerically, we have checked the validity of this conjecture for all odd up to 101. 

Proof of Theorem 12. 5[ The proof is similar to the above proof of Theorem l2.3l without 
using any information on the sign of Co,2Ar-i- D 

Conjecture 12.61 then follows from Conjecture 14.91 by including the information on the 
sign of Co,2Ar-i in the proof. 
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• y = g_ (x) 



Figure 9. Images of the boundary of the square [—1, 1] x [—1, 1] under the map H, shown 
(a) for 7 = 1/4, and (b) for 7 = 0.258 . . . The "vertical" strips V_, Vq and V+ are mapped 
on "horizontal" strips H{V-), H{Vo) and iJ(V+). Their intersections define the elements of 
the first-level partition, for instance /_._ —V-H HiV-) is the lower-left squarish domain. 
Intersections of the light curves, obtained by iterating the map twice, define the elements 



of the second-level partition. 



A Small Coupling and Symbolic Dynamics 

In this appendix we sketch the proof of Proposition 12.11 on the continuation of equilibrium 
points from the uncoupled limit. The equation -|- ^(x^+i — 2x„ -|- = satisfied 

by the stationary points of the potential can be rewritten as {xn+i-,yn+i) = H{xn,yn) 
where H is the map of the plane defined by 

H{x,y) = (2x-^fix)-y,x) . (A.l) 

This map is invertible and its inverse can be simply obtained as = S o H o S where 
S{x,y) = {y,x). In addition, we have H{—x,—y) = —H{x,y). 

We proceed to a similar construction as in [Kee87j to show that, when 7 ^ 1/4, the 
map H has an invariant set on which it is conjugated to the full shift with 3 symbols. This 
implies that the map admits at least 3^ points of (not necessarily minimal) period N for 
any N. Since the potential is a polynomial of degree 4 in variables, it can have at 
most 3^ stationary points, and thus we will have obtained the exact number of stationary 
points. 

We start by defining a collection of strips in the square [—1, 1]^ with good dynamical 
properties. The map x ^ g{x) :=2x — 2'y~^f{x) + 1 leaves —1 invariant and is strictly 
increasing on [—1, —zq] where zq '■= ^J{l — j)/3. If in addition, the condition g{—ZQ) ^ 1 
holds, then g is a bijection on [—1, gZ (1)] where gZ^ is the inverse of g with range [—1, —zq] 
and we have 

h(^{{x,-1): gZ\-l)^x^gZ\l)}) ={{x,gZ\x)): - 1 ^ x ^ l} . (A.2) 
Similarly, provided that the condition g{—ZQ) ^ 3 is satisfied, we have 

h(^{{x,1): gZHl)^xi^gZ\S)}^ ={{x,gZHx + 2)): - 1 ^ x ^ l} . (A.3) 
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The condition g{—ZQ) ^ 3 is indeed equivalent to 7 ^ 1/4. By monotonicity, it results 
that under the condition 7^1/4 the "vertical" strip 

V-:={{x,y): gZHy)^x^g-Hy + '^), -1 ^ y ^ l} , (A.4) 

is mapped onto the "horizontal" strip S{V-). By symmetry we also have that if(V+) = 
5(V+) where V+ = —V- is the opposite vertical strip. 

Furthermore, g is strictly decreasing on [— 2:o,zo] and the conditions g{—zo) ^ 1 and 
gi^o) ^ ~1 imply that this map is a bijection on [g^^ {1) , g^^ {—!)] where g^^ is the inverse 
of g with range [zq, —zq]. The inverse g^^ has the symmetry g^^{2 — x) = —g^^{x) and 
the condition g{zQ) ^ — 1 is equivalent to g(— zq) ^ 3. Similarly as before this implies that 
under the condition 7 ^ 1/4, the symmetric vertical strip 

Vo:={{x,y): g^\y + 2) ^ X i;: g^\y), -1 ^ y ^ l} , (A.5) 

is mapped onto the horizontal strip 5(Vo). 

With these strips provided, the rest of the construction is standard. Given any finite 
word LJ_„ • • • uJn+i G { — ,0, (^n ^ Q) consider the set 

n+l 

k=—n 

These sets are all non-empty by the properties of the strips above. Then, for any infinite 
sequence {oJn} S {— )0, +}^, the intersection Hfcgz H^^l^uk) is ^^Iso nonempty. Since the 
sets /a;_„ --aj„+i are pairwise disjoint, we conclude that the map H has at least (and thus 
exactly) 3^ points for which the forward and backward orbit is contained in [—1,1]^. In 
particular it has exactly 3^ points of period N for any N. This shows that ^*{N) ^1/4 
for any N ^ 2. 

The condition g{—ZQ) ^ 3 which ensures that the strips have the desired properties is 
not optimal and can be improved. To see this, we consider the strips V^^, instead of Vuj^. 
in ()A.6jl where V!_ is defined by 

V'_:={ix,y): - 1 ^ x ^ -zo, g{x) - 2 ^ y ^ g{x)} , (A.7) 

where = —V'_ and where 

V^:={(x,y): - zq ^ x ^ zq, g{x) - 2 ^ y ^ g{x)} , (A.8) 

Note that we have V[j^ D Vuj^ ■ 

Similarly as before, a sufficient condition for the existence of 3^ points with bounded 
orbit under H is that the nine basic sets I^jq.ui are pairwise disjoint. By symmetries, 
the latter is equivalent to the condition that the right boundary of crosses the upper 
horizontal strip S{V'^). A simple analysis shows that this is equivalent to the inequality 
(weaker than g(— zq) ^ 3) 

g{-zo)-2^gZH2-zo) (A.9) 

which holds provided 7 ^ 0.258 . . . , i.e., we have 7*(A^) ^ 0.258 ... for any N ^ 2. 

As 7 increases, some symbolic sequences will be pruned, resulting in saddle-node bi- 
furcations of stationary points. The first configurations to disappear are asymmetric ones, 
such as (1,1,..., 1, —1) and (1,1,..., 1, 0) and their images under the symmetry group G. 
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Orbit label 


Fixed-point set 


Reduced dynamics 


I± 


(x, X, X, x) 


X = X — x^ 


A 


{x, X, —x, —x) 


X = (1 — 7)x — x^ 


B 


{x, 0, —X, 0) 


X = (1 — 7)x — x^ 


^(2) 


(x, —X, X, — x) 


X = (1 — 27)x — x^ 


Aa 


{x,x,y,y) 


X = (1 - ^7)x + ^72/ - x^ 






y = \-ix + (1 - \-i)y - y^ 




{x,y,x,y) 


X = (1 — 7)x + 7y — x^ 






y = 7X + (1 -7)y - y^ 




{x,-x,y,-y) 


X = (1 — |7)x — ^72/ — x^ 






ij = -\-ix + (1 - \i)y - y^ 


da, db,. . . 


{x,y,x,z) 


X = (1 — 7)x + \^y + — x^ 






2/ = 7X + (1 - 7)2/ - y^ 






i = 7X + (1 — 7)z — 



Table 6. Orbits, corresponding fixed-point sets and reduced dynamics in the case TV = 4. 

B The case = 4 

In this appendix, we briefly describe how we obtained the bifurcation diagram of Figure [5] 
for the case = 4. 

We start by determining the isotropy groups of the stationary points present in the 
uncoupled case 7 = 0, and their fixed-point sets. They are of 10 different types, 8 of which 
we show in TablelH The two types we do not show are the fixed-point set {0, 0, 0, 0}, which 
obviously contains only the origin, and the orbits labelled Aaa, which have no symmetry, 
and thus a fixed-point set equal to X. 

Table [6] also shows the form taken by the equation x = — VV^(x) when restricted 
to the fixed-point set. These equations are then analysed in order to determine their 
equilibrium points. The analysis of the dynamics in one-dimensional fixed-point sets is 
straightforward, while the two-dimensional equations require a bit more work, but do not 
present particular difficulties. Let us briefly describe the case of the three-dimensional 
fixed-point set (x, ?/, x, z). In rotated variables u = {y + z) /2, v = [y — z)/2, the equations 
become 

X = (1 — 7)x -I- 7ii — x^ , 

tt = 7x + (1 — 7)ti — u{u^ + 3d^) , (B.l) 
i) = {1 — j)v — v{3u^ + v^) . 

Equilibrium points with v = Q correspond to the fixed-point set {x,y,x,y), and have 
already been analysed. We thus assume v ^ 0, which implies Su^ -|- f = 1 — 7 and allows 
to eliminate v. We are left with two equations for the stationary points, namely 

(1 - 7)x 7M - x^ = , 
7X - 2(1 - 7)n Sn^ = . 
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If j; = 0, then necessarily u = and we are left with the already studied points of the form 
(0, y, 0, — y). We thus assume x ^ and eliminate u from the system. It is convenient to 
introduce the variable w = 2(1 — 7 — x^) instead of x. Then u = —wx/2j, and w has to 
satisfy the fourth-order equation 

- 2(1 - 7)^3 + 27^(1 - j)w + 2-f^ = . (B.3) 

The condition ^ yields the additional requirement — 2(1 — 7)7/;^ + §7^(1 — 7) ^ 0. 
Together with (|B.3p . this can be seen to be equivalent to the condition 

37' . ^ 

— ^ 1 - 7 • B.4 

w 

Numerically, one observes the following properties: 

• If < 7 < 0.2684 . . . , Equation ()B.3P has four distinct real roots, two of them being 
positive and two negative. 

• For 0.2684 ... < 7 < 0.4004 . . . , Equation ()B.3P has two real roots, which are both 
positive. 

The negative roots always fulfil Condition ()B.4p . Examining ()B.3P perturbatively in 7, 
one finds that they are of the form w = —7 + 0(7^) and w = —7^ + 0(7'^), and correspond 
to the branches labelled da and db in the bifurcation diagram of Figure [5l 

The positive roots, on the other hand, do not always fulfil Condition (IB.4p . In fact, 
one solution exists for < 7 < (3-v/2 — 2)/7, while the other one bifurcates with the 
^(^^-branch for 7 = 2/5. A perturbative analysis for small 7 shows that they converge, 
respectively, to the points (1, 0, 1, —1) and (0, 1, 0, 0) as 7 ^ 0. 

Finally, one can also compute the determinant of the Hessian of the potential around 
the stationary points of the form {x,y,x,z). This determinant can be put into the form 

^ [3w - 4(1 - 7)] [(1 - 7)^1- - 37^] [3(1 - 7)w;2 ^ _ 2^)^^ + g^a^i _ ^)] . (b.5) 

The roots of this expression correspond to bifurcation points, and can be combined 
with ()B.3P to check that all bifurcation points have been found (and yields more pre- 
cise estimates of the bifurcation values). This expression can also be used to check the 
index of the stationary points. 

This analysis is not completely rigorous because we have not analysed in detail the 
branch labelled Aaa, which has no symmetry at all. The bifurcation diagram is based on 
the fact that the only bifurcation of already known branches producing stationary points 
without symmetry is the bifurcation of the Aa-branch at 7 = 1/3. A local analysis shows 
that this is indeed a pitchfork bifurcation, producing the right number of new stationary 
points. 
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